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Abstract 

We derive the matching conditions for cosmological perturbations in a Fried- 
mann Universe where the equation of state undergoes a sharp jump, for in- 
stance as a result of a phase transition. The physics of the transition which 
is needed to follow the fate of the perturbations is clarified. We dissipate 
misleading statements made recently in the literature |ij concerning the pre- 
dictions of the primordial fluctuations from inflation and confirm standard 
results. Applications to string cosmology are considered. 
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I. INTRODUCTION 



The evolution of metric perturbations in a Friedmann Universe is well known (see, e.g. 
and refs therein). They can be classified into scalar perturbations, which couple to the 
matter density inhomogeneities, and freely propagating gravitational waves (we shall ignore 
here the vectorial perturbations). In inflationary models [3] the perturbations of both types 
are produced from primordial quantum fluctuations, and their spectra after inflation are 
nearly scale invariant [4,5]. However the amplitudes of the gravitational waves and scalar 
perturbations are significantly different at all scales after the transition from the inflationary 
to the matter dominated stage. 

Let us, for instance consider a Universe which evolves from a phase (-) with the scale 
factor a_ oc t p ~ to a phase (+) such that a + oc t p+ . If the transition from the (-) and to 
the (+) phase happens at time t t , then the ratio of the amplitude h + of the long wavelength 
gravity waves some time after the transition (t ^> t t ) and directly before it h~ (at t = t t — 0) 
is: 

f^om (i.D 

while for the gravitational potential $, which characterizes scalar metric perturbations, the 
result is 

£ - . (1.2) 

We see that for a strongly inflating (-) phase (p_ ^> 1) the amplification of density pertur- 
bations can be very big whereas gravity waves are not amplified. 
This standard result has however been challenged recently in 

We think that there is nothing strange about this difference in the amplification of the 
gravity waves and density perturbations. Of course, eq. ( |1.2|) doesn't imply that the gravita- 
tional potential jumps at the moment of transition. What happens is that the gravitational 
potential $ which can be small immediately before and after the transition is distributed at 
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t = tt + O into two modes, both very large in amplitude, one which decays, the other yielding 
( |1 . 2| ) (the derivative <3? of the gravitational potential on the other hand can strongly jump 
at t t in distinction from the gravity waves for which h is continuous). As we will show the 
result depends on the physics of the transition. 

In order to make this statement quantitative we derive the matching conditions for 
cosmological perturbations in models with sharp transitions. The results will be applied to 
the simplified model of inflation considered in |]J which uses a sharp jump of the equation 
of state at the end of inflation instead of evolving it rapidly but smoothly as in e.g. (see 
also e.g. H0). We think that the results obtained in |] imply the strongly scale dependent 
spectrum for scalar perturbation in the relevant for COBE scales. We do not agree with the 
author of paper ]IJ who concluded from his consideration that the spectrum is flat taking the 
parameter 71 to be equal 2 on p. 7166, while it should be taken 71 <C 1- It will be shown in 
Sect. |V B| that the root for the misleading result lies in inappropriate matching conditions 



which have no physical justification. We on the other hand shall make the matching on 
hypersurfaces of constant energy for physical reasons and clarify and confirm the standard 
results O-O . 

The considerations in this paper actually go beyond clarifying this point. They can be 
applied to study the evolution of cosmological perturbations in models where the transition 
from one stage of evolution to an other is not completely specified. We clarify the physical 
features of the phase transition which are of importance for the evolution of cosmological 
perturbations. As a particular example we consider the fate of the growing mode in string 
cosmology after transition from superinflation to the Friedmann era (see Sect. |V|). 

The paper is organized as follows. In Section |J we recall the standard results. In 
Section [TTT] the general matching conditions in models with sharp transitions are derived. 
We also discuss there the choice of the hypersurface of transition. In Section [IV] we write 
the matching conditions in two particular coordinate systems: the longitudinal and the 
synchronous reference frame used in UJ. Finally, in Section |VJ we apply our results to 



concrete examples. Conclusions are drawn in Section [VI. 



II. BACKGROUND AND PERTURBATIONS 



Consider a quasi-isotropic spatially flat Friedmann Universe with metric 

ds 2 = a 2 {(l + 2(p)dr] 2 - 2B ti da?dri - [(1 - 2ip)5 lk + 2E )ik + h ik \dx l dx k } (2.1) 

where a ir\) is the scale factor, a comma means the derivative with respect to spatial coordi- 
nates and the four functions 0, ip, B and E characterize the metric perturbations of scalar 
type. The tensor hi k , which is taken transverse and traceless {h lk k = h\ = 0) corresponds to 
gravity waves (all spatial indices will be raised with S lk ). We ignore here the vector pertur- 
bations. Under coordinate transformations rj — > fj = r\ + £°, x % — > x l = x % + £> l , the scalar 
metric perturbations transform as: 

0^0 = 0-^°-^°' ; B^B = B + e-e'; 

^-^ = -0-^° ; E^E = E-i- (2.2) 

where Ti = a' /a = Ha and a prime denotes a derivative with respect to conformal time rj. 

From ( |2.2j ) we see that two simple gauge invariant quantities, characterizing the scalar 
perturbations, can be built, e.g. 

$ = (j) + -[{B - E')a]' , * = i[> - H(B - E') . (2.3) 

(X 

In the longitudinal gauge (Be — Eg — 0) these gauge invariant variables coincide with the 
metric perturbations. Therefore the calculations in this coordinate system are completely 
the same as those in terms of the gauge invariant variables $ and ty. As for the gravity 
waves they are gauge invariant: hij = hij. 

If we consider the perturbations with scales much bigger than the Hubble radius H -1 , 
then the solution of Einstein equations in the long wavelength approximation can be written 
down as (see e.g. 0): 

h = D g + S g J^ (2.4) 

for gravity waves and 
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$ = V = S a — + D 8 (-[adt) (2.5) 
a \a J J 

for scalar perturbations in the case of matter an adiabatic perfect fluid or a scalar field. 
Here a dot denotes a derivative with respect to cosmic time t = J adr], and H = a/ a. The 
constants of integration D, S are fixed by physical considerations, e.g. by imposing that 
the perturbations be in their quantum ground state when they are inside the Hubble radius 
during an inflationary era. We note that the solution ( |2.5| ) can be obtained if we integrate 
the equation for the quantity (see, for instance PJ5jflp!0|1) 

which is "conserved" in the long wavelength approximation. 

Let us pause onto the model described in the Introduction, where the scale factor evolves 
as a oc t p ~ before the transition and as a oc t p+ after it. Then using eqs. fl2.4|) , fl2.5|) one gets 
that h~ — > D g , $~ — > D s /(1 +P-) and after the transition after the subdominant mode has 
decayed: h + —>■ D g , $ + — >• D s /(1 + p + ). From here immediately follow the results ( [L.1D , 
( |1 . 2\ ) . Since during inflation 1 + p_ oc H 2 /H > 1 we see that the scalar perturbations can 
be very strongly amplified in the course of the transition. 

The standard derivation described above was criticized in |TT| on the grounds that the 
long wavelength approximation used breaks down for scalar perturbations of cosmologically 
relevant scales at the moment of the transition. It is true that the condition of applicability 
of the long wavelength approximation for scalar perturbations can be significantly different 
from simply imposing that the wavelength be bigger than the Hubble scale (see || and 0). 
Therefore we believe that it is quite useful to clarify how the perturbations evolve in models 
with a sharp change of the equation of state and discuss the features of the physics of the 
transition which are relevant for cosmological perturbations. 
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III. MATCHING CONDITIONS 



Suppose that the stress-energy tensor, which governs the evolution of the metric ( |2.1|) 
via Einstein's equations undergoes a finite discontinuity on a spacelike hypersurface £ de- 
termined by the equation q{r],x l ) = const., where q is 4-scalar. In a slightly inhomogeneous 
Universe it can be decomposed into a "homogeneous" part qo{if) plus a small inhomogeneity 
5q{rj,x l ), that is q = qo + 5q. Under coordinate transformations 77 = fj — 77 + £° etc., 5q 
changes as 

Sq^5q = 5q-q' o e ■ (3.1) 

We shall consider below specific examples for the function q. The hypersurface £ divides 
the manifold into two regions, (— ) and (+) which should be matched in such a way that 
the induced 3-metric on E and its extrinsic curvature be continuous (in the absence of 
surface layers) [f[^| . The simplest way to do the calculation is to go to the "tilde coordinate 
system": rj — > fj, x % — > x 1 , where the equation for £ becomes fj = const.. It means that in 
this coordinate system 5q = and eq. ( |3.1|) determines the transformation: 

£° = 4" > £ - arbitrary . (3.2) 
% 

From this, one first immediately gets that for the background and gravity waves the matching 
conditions impose that the scale factor a, and their first time derivatives be continuous 
on £ in any coordinate system. For the scalar perturbations we can easily calculate the 
extrinsic curvature 5Kj in the "tilde coordinate system" and the matching conditions for 
them read: 

[$± = [E]± = , (3.3a) 
[SK]]± = -~[8}{H4> + -00 + (B- &)'% = , (3.3b) 

(X 

where [if)]± = ip + — ip- etc. From eq. ( |3.3b| ) it follows that (H(fi + tj)') and (B — E') should 
be separatly continuous. Expressing if}, B etc... in terms of if), B etc... in the original 
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coordinate system with the help of (|2.2|) where £° should be taken from (|3.2|) , we finally get 
the matching conditions on E: q Q + 5q = const, in an arbitrary coordinate system: 

[il> + H6q/<&]± = 0, (3.4a) 

[E-£] ± = 0, (3.4b) 

[B-E' + 6q/q' ] ± = 0, (3.4c) 

[H(P + ij' +(H' -H 2 )5q/q' }± = . (3.4d) 



Eq. (|3.4b ) is empty since £ can always be chosen so that it is satisfied. 



Now it is time to specify E, that is the scalar q. Suppose matter is an adiabatic perfect 
fluid with stress-energy tensor = (p+e)u^u u —p5^, and that the equation of state p = p(e) 
undergoes a sudden change at some moment of time. Then it is clear that on E the energy 
density e should be constant (that is, E: q = e + Se = const.), since the pressure depends 
only on e. The background energy density 6q and its perturbations 5e can be expressed in 
terms of the scale factor a, and the metric perturbations <fr, ip etc... via the linearized — i 
Einstein eqs. and thus the eqs. ( |3.4j ) can be written down entirely in terms of the metric 
perturbations. We will do it in next Section for various gauges. 

If matter is a scalar field ip with stress-energy tensor = ip'^(p )U — 6^(^(p' x (p^ — V(ip)) 
which suddenly decays into ultrarelativistic particles then E should be taken to be the 
hypersurface e = \i{)' x ip t \ + V((p) = const., where e is the "effective" energy of a scalar field 
considered as a "perfect fluid" with 4-velocitie = ip tft / ' (ip' x (p^\) l l 2 . 

IV. LONGITUDINAL VS. SYNCHRONOUS GAUGE 

A. Longitudinal gauge 

The longitudinal gauge conditions: Et = Be = 0, define a unique coordinate system. In 
this coordinate system the matching conditions (|2.4j) on the hypersurface E : e + $ £ e 
const simplify to: 
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hrl± = o, 



[Ufa + i// t + H'5e e /e ]± = 



(4.1) 
(4.2) 



Now we will use the background — Einstein equation to express Eq and e' in terms of the 
scale factor and its derivatives: 



m 2 



'0 2 



(4.3) 



where k = 8nG, and the i — k (i ^ k) (if is supposed that 5TJ: = for i ^ k together with 
— linearized Einstein equations in the longitudinal gauge (see Sec. in 0) 



ipi = <j>t, bet 
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KCL Z 



(4.4) 



Then substituting these eqs. in ( |4.2| ) we arrive at the following two independent matching 
conditions (the third one becomes redundant): 



0. 



U 3H'-H 2 







J ± 



(4.5a) 
(4.5b) 

). The second 



where Q is the "conserved" quantity defined by (|2]6|) (recall that $ = 
term proportional to A(p£ in fl4.5b| ) can be neglected for the perturbations with scales bigger 
than the Hubble scale if -1 . Then we see that for these long wavelength perturbations the 
matching condition ( |4.5 b| ) reduces to the widely used conservation law for Q. It is worth to 
mention that for long wavelength perturbations the continuity of <pi and Q also implies that 
the 3-velocity u\ should be continous on E. 

We conclude this paragraph by noting that all Eqs. (|3.1|)-(|3T4]) can be turned into 
relations between gauge invariant quantities by changing ip£, into <fii into $ etc. This 
remark allows to readily translate these Eqs. in any other reference system. 



B. Synchronous gauge. 



The synchronous frames are defined by imposing the gauge conditions: (J) s = B s = 0. It 
is well known that these conditions do not fix completely the coordinate system. Indeed, 
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under the coordinate transformations 

Vs^Vs = V + ~, x* -> x* = x { + A'' /— (4.6) 

a J a 

where A is an arbitrary function of the spatial coordinates, the gauge remains synchronous 

(0 S = B s = 0) and (see eqs. Q). 

</> s -> & = i>s + ; E s —> E g = E g — \ — . (4.7) 

a J a 



The function A and the constant of integration in (|4.7| ) correspond to the so-called fictitious 
modes which do not lead to any physical inhomogeneities. The contraint eqs. Q4.4|) in 
synchronous frames go into (see, e.g., @): 

"I 



HQ? 



3 A^ s +nE' s )-w s 



(4.8) 



with the help of eqs. ( |4.3| ) , ( ^4.8| ) we can easily deduce from ( |3.4| ) the independent matching 
conditions in synchronous coordinate systems on the hypersurface £ : £q + 5e s = const, 



. (4.9) 



fok + HE' S ] ± = 0, [HE' S - n , l _ n2 [l^s + HE' S ) - W s 

First we note that the particular combinations of the metric variables entering ( f4.9p do not 
depend on the fictitious modes and we can work directly with ( f4.9| ) to find the amplitudes 
of the physical modes without specifying a particular synchronous frame. Second, if we 
substitute the expressions for ip£, En in terms of <fie in (|4.9|) we immediatly arrive to ( (4.5|) . 
So the matching conditions (fO[ ) are completely equivalent to (|4.5|) as it should be. 

We see from (^4.9| ) that, contrarily to what is claimed in |Q (begining of |V|, p. 7164), the 



matching conditions in a general synchronous frame can not be reduced in general to the 
continuity of metric and its first time derivatives. The reason beeing that the hypersurface 
E that we choose on physical grounds as the surface Eq + 5e s = const is not necessarily the 
surface r] = const. However, using the remaining coordinate freedom ( |4.6| ) we can pick up 
the particular synchronous frame in which the hypersurface £ : Eq + 5e s = const, coincide 
with the hypersurface r] s = const. It is easy to see from (|4.8|), ( f4.9| ) that this is achieved if 
we take 5e s \y, = 0, or in explicit form 
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. (4.10) 



Then in the particular synchronous frame defined by ( |4.10| ) the matching conditions ( |-1.9| ) 
reduce to the continuity of if) s , if)' s , E' s . The quantity E s itself can always be made continuous 
on E since there is an extra constant of integration in synchronous frames. These continuity 
conditions together with fl4.10|) completely fix all the constants of the problem. For the long 
wavelength perturbations when the term A(^ s + 7iE' s ) can be ignored, eq. fl4.10|) simplifies 
to t/^Ie = 0. Since the perturbations of the 3- velocity 5u l s are proportional to (it follows 
from — i Einstein eqs), we conclude that for long wavelength perturbations the coordinate 
system which satisfies ( |1.10| ) is the particular synchronous frame which is comoving on 
the hypersurface E, that is 5u\ |s = 0. Finally, we would like to stress that the matching 
conditions discussed here for longwave perturbations automatically imply the conservation 
of ( during the transition. 



V. EXAMPLES 

A. From a pure de Sitter to a dust dominated era 

To start with we would like to consider some artificial example to demonstrate that we 
can generate physical inhomogeneities as a result a of transition even if we start from a pure 
de Sitter Universe where it is known that the scalar metric perturbation are exactly equal 
to zero. As we will show the result is determined by the physics of the transition. First we 
recall some useful formulae which relate the metric perturbations in a synhronous frame to 
the gauge invariant gravitational potential which can be easily obtained by integrating ( p.3|) 
if we use there the conditions B = (f> = and take into account the fact that \P = $: 

ip s = $ + — J ®adr], E s = - J — J Qadr] . (5.1) 

In a pure de Sitter Universe $ = 0. The constant of integration A in ( |5.1| ) (for instance, 
?p s = A—) corresponds to fictitious perturbations. If we fix it to be Ao 7^ then we pick 
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up (up to trivial time independent transformations) a particular synchronous coordinate 
system from the wide class of synchronous flames which are related one to each other by the 



transformations ( |4.6| ). Let us assume (without any physical justification) that the transition 
from the de Sitter to the dust era happens along the hypersurface S : 770 = const, in this 
particular reference frame characterized by Aq. The solution for the gravitational potential 



in a dust dominated universe (a m oc rj m ) is (see( |2.5D ) 

$ = An + ^ , (5.2) 
Vm 

where both constants of integration D m and S m correspond to physical (nonfictitious) modes. 
Then using ( |5.1|) one can easily derive the solutions for the metric perturbations in the 
synchronous frame during the de Sitter stage: 

^ = X H d . s , K = ~> ( 5 - 3 ) 

where the Hubble constant H^.s = — does not depend on time in a pure de Sitter Universe, 
and correspondingly in the dust dominated era: 

'P,= lD m + 2 Jf,E', = -^r lm + 1 -^-^, (5.4) 

where F m is the constant of integration corresponding to the fictitious mode. From the 
continuity conditions for x/j s , ip' s and E' s which imply the continuity of the metric and extrinsic 
curvature if E : 77 = const, (the continuity of E s fixes only the extra constant of integration 
in ( |5.2j) ) one immediatly gets that 



An = -^oHd.s , [~t] = ~T^oHd.s , F m = . (5.5) 



Since T> m corresponds to physical inhomogeneities we see that as a result of the transition 
from a pure de Sitter to a matter dominated Universe one produces real scalar perturbations. 
Of course immediatly after de Sitter stage the gravitational potential is still zero: $e+o — 0, 
but it consists of two modes, one which decays, while the other survives. The final result 
for the amplitude of the physical perturbations is determined by the hypersurface of decay 
of the de Sitter Universe (by Ao in example considered here). 
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Now, in a pure de Sitter space there is nothing which could help us to fix the hypersurface 
of transition to a matter dominated era. For instance, any hypersurface is an hypersurface 
of constant energy. So the physical result is uncertain. Therefore the small deviations from a 
pure de Sitter stage play a crucial role in predicting the resulting inhomogeneities, since they 
actually define the physics of the transition. The most natural possibility (in the absence 
of further concrete model) is to take as £ the hypersurface of constant energy density as we 



Now we consider the model studied in |]]. Namely, we assume that the Universe went 
through three stages: at first there was inflation (i) with scale factor oc (— r^)" 1 " 7 , where 
< 7 < 1, then a radiation dominated era (r) with a r oc r] r and finally a dust matter 
dominated stage (m): a m oc rj m . Let us denote the times of transition from i to r stage 
as r)i and correspondingly from r to m stage as 1]2 in terms of conformal time r) r in the 
r-stage. Under the above parametrization of the scale factors, conformal time rj jumps 
during the transitions and to relate, for instance rji to r] r at this point we can use the 
continuity condition for Ti. The time r\ can be easily made to be continuous by shifting r] r 
and rj m by constant factors. However it has no impact on the final results and to simplify 
the formulae we prefer to use the parametrization above. Also from the very beginning 
we restrict ourselves only to inhomogeneities with scales bigger than the Hubble radius at 
rj = 772, that is (^772) <C 1, where k- is the comoving wavenumber of planewave perturbations. 
Then for such perturbations we can write the solutions for the gravitational potential $ in 
the long wavelength approximation during the inflation and radiation dominated stages as: 



did. 



B. Prom inflation to a radiation and then dust dominated era 



(5.6a) 




(5.6b) 
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We used Planck units and Hf is the value of Hubble constant at the moment of horizon 
crossing by the perturbation with the wavenumber k during inflation. Imposing the pertur- 
bation to be in their quantum grond state during inflation fixes the amplitude in ( [5.6a ) Q 



(there is no disagreement in the literature on this point) and we ignored the decaying mode 
which will be completely irrelevant at the end of inflation. In the radiation stage we keep 
in (|5.6b|) the first correction ~ (cor]) 2 , where to = k/\/3, to the constant nondecaying mode 



since at leading order this mode has a vanishing derivative. 



The solution ( |5.6b| ) can be obtained from well known exact solution 0JI| in the r-stage, if 



we expand it in powers of wrj. Using (|5.6| ) one easily gets with the help of (|5.1|) the following 
expression for ip s and E' s in a synchronous frame: 

^ = ^^#f + (l + 7M-^) 7 , 

^' = -^f(-^)-^(-^) 1+7 , (5-7) 
in the inflationary stage (i) and 



| (l - >,,)*) Dr + I 



in the r-stage, where the constants of integration Fi, F r correspond to fictitious modes. 

The potentials in a dust matter dominated Universe were given before (see formulae ( |5.2| ) . 
( |5.4|) ). Now if we claim that ip, ip f , E' are continuous two on the transition hypersurfaces, 
r] r = r/i and r\ r = r] 2 , then we get 6 eqs. for seven unknown constants D r , S r , D m , S m , F iy 
F r , F m . Clearly we need one extra condition. The one imposed in [J] is that in the dust 
dominated era the synchronous coordinate system be also the comoving one. This fixes the 
constant F m to be zero. Then a straightforward calculation gives us the following result for 



the constant D m : 



m ■ 

3(l + 2 7 ) (, , 2 + T/ _._ , 2 , . 



D m = + -^WiWiirj (5-9) 



which agrees with [0]. 

13 



Recall that D m is the gravitational potential in the m-stage. Since 7 < 1 we see from 
Q5.9D that for cosmologically relevant scales for which ^(tor]2) 2 (r]2/r]i) 2 1, the spectrum 
Q5.9| ) is by no means scale invariant (for which D m doesn't depend on k). We do not agree 
with the author of paper |1| who concluded on the base of consideration similar to the above 
one that it implies a scale invariant spectrum in relevant for COBE scales. The reason for 
disagreement is that he took, as we think mistakenly, in formula for T> on the page 7166 
IjH the parameter 71 to be equal 2 (it corresponds to 7 = — 2 in our notations) instead of 
71 <C 1. Since the spectrum ( |5.9|) doesn't resemble at all the standard result (see 0) it is 
quite interesting to find out the reason for this desagreement. 



As we argued in Sect. |T|, the physics of the transition tells us that it should happen 
on the hypersurface of constant energy. This implies that, for instance, in an arbitrary 
synchronous coordinate system the matching conditions should be (^4.9| ). They are equivalent 
to the continuity conditions for ip s , ip' s , E' s only in the particular synchronous frame which 
for long wavelength perturbations is close to the comoving one at the moment of transition. 

The precise condition which fixed this coordinate system and, accordingly, the amplitudes 
of the fictitious modes F is given by ( |4. 1 0| ) . If we decide the synchronous frame to be 
comoving one in the m-stage (F m = 0) then this coordinate system satisfies the requirement 
( |4.10| ) with very good accuracy for kr]2 C 1 on the hypersurface of transition from radiation 
to matter dominated era. The function F r which specifies this particular synchronous frame 
II in the r-stage is u F r ~ u 2 (7]2) A D r . However, the frame II is very far from the comoving 
one on the hypersurface of the first transition from the i to the r stage. Actually in the frame 
I which is comoving on the hypersurface of the first transition the amplitude of the fictitious 
mode is l F ~ u 2 (r)i) A D r . Since ll F/ T F ~ 772/171 3> 1, we see that the frames I and II are 
very different. As one can check this difference is quite important for the perturbations with 
0^72 ) 2 (^/^i) 2 3> 1. How then should we proceed to get the correct answer for the final 
amplitude of scalar perturbations ? The simplest way is just to use the matching conditions 
( |4.9[ ) which do not involve the fictitious modes oc F and permit us immediately to write 
4 eqs. for 4 unkonwn constants D r , S r , D m , S m . However, if we want to insist on the 
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continuity conditions for ip s , ip' s and E' s then we can do it but, after the matching in the 
particular frame I at rj r — r\\ we should go to the other synchronous frame II to make the 
matching at r\ r — r]2- This procedure fixes all the constants of the problem. As one can easily 
check both ways lead to the standard result for the final amplitude of scalar perturbations: 

3 

D m = -^, (5.10) 

for 7< 1, which is very different from ( |5.9|) , while for longwave gravity waves h ~ H\ . 

Finally there is noting strange about the fact that the final amplitude of scalar pertur- 
bations can go to infinity when 7 — > despite the fact that perturbations during inflation go 
to zero (see ( |5.6a|) ). As it was shown in Sect. |V A| as a result of the decay of a pure de Sitter 
stage we can get whatever we want unless the physics of transition is specified. In the model 
considered above the hypersurface of transition becomes more "inhomogeneous" (we mean 
its extrinsic curvature) when we approach a pure de Sitter space. The formulae ( |5.10|) has 
however a limited range of validity since it was obtained in the linear perturbation theory. 
Thus we can believe it only if 7 ^> Hj ($ ~ D < 1). We hope that the above consideration 
clearly demonstrates that the calculations which are very straightforward in terms of gauge 
invariant variables (or in the longitudinal gauge) can become sophisticated in synchronous 
frames. 



C. Fate of the "decaying" mode (applications to string cosmology) 

Finally we want to consider some consequences of the obtained results for string cosmol- 
ogy. It is well established (see for instance, |13[) that in the models with a dilatonic field 
in the absence of nonperturbative potential there is a stage in the evolution of the Universe 
which is described by "superinflationary" expansion in the "Brans-Dicke frame". It is usu- 



ally assumed that this stage should serve the same purposes as ordinary inflation |13[]. In 
the "Einstein frame" this superinflationary stage corresponds to a contracting Universe and 
the "decaying" mode in Q2.5|), 
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$ = S— , (5.11) 
a 



(we consider only long wavelength inhomogeneities) grows in this contracting Universe |L4 
There is therefore the danger that this mode will eventually turn in to a nondecaying mode 
and create big inhomogeneities in the final Friedmann Universe thus completly invalidating 



the original model. It was found in [|14]] that the mode ( |5.11| ) has peculiar properties which 



give us the hope to avoid this danger. However, the question of what will happen with it 
after the transition from the superinflationary to the Friedmann stage was not clarified. The 
problem is that there are no realistic models which solve the "graceful exit" problem and 
allow us to go smoothly from the superinflationary to Friedmann era fl5| . In the absence 



of such a model we will assume that the transition happens as a result of a change in the 
effective equation of state along a hypersurface of constant energy. On the basis of the results 
obtained above it turns out to be enough to predict the fate of the mode (|5.11|) (which is 
growing in superinflationary Universe) after transition to Friedmann era. 

As we have seen before the nondecaying mode of the gravitational potential in an ex- 
panding Universe after the change in the equation of state is redistributed in comparable 
proportions between nondecaying and decaying modes. If it would be also the case for the 
decaying mode the situation would be quite unfavorable for the superinflationary models. 
However substituting ( [5.11| ) in ( |4.5| ) we immediately find that 



S + = S- (5.12) 

to leading order in (fc^)-expansion. A more careful analysis shows that the nondecaying 
mode will be generated only at the next order in (kr]), that is D + ~ (kr]) 2 S~. Thus we 
see that the growing mode of scalar perturbations in superinflationary Universe will be 
practically entirely converted in the decaying mode in Friedmann Universe. This result is 
general and doesn't depend on the concrete model for transition from superinflationary to 
Friedmann era. The only restriction is that the transition happens on a hypersurface of 
constant energy. 
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VI. DISCUSSION 



We derived the matching conditions for cosmological perturbations during a sharp change 
in the equation of state in gauge invariant form and in different gauges (including the 
synchronous one) assuming that the transition from one stage to the other happens on the 
hypersurface of constant energy. This last restriction is quite natural and it is dictated by the 
physics of the transition. Using it we confirmed and clarified the standard result concerning 
the ratio of density perturbation to gravity waves predicted by inflationary models. 

We showed that the matching conditions used in OJ do not satisfy this requirement and 
therefore lead to results for the spectrum of density perturbations at odds with the standard 
results. 

The results obtained go beyond the model usually considered and can be applied also 
study to the fate of cosmological perturbations in models where not all of the details of the 
physics of the transition are specified. The only important piece of physics which we used is 
that the transition happens on a hypersurface of constant energy. The results were applied 
to follow the fate of the growing mode in superinflationary cosmology . 

A final remark which we would like to make concerns the "best" gauge for cosmological 
perturbations. As far as we are aware the question about "what is the best" gauge or 
the advantage of the gauge invariant formalism is still under debate in the literature. Our 
(subjective) point of view is quite simple. There are of course no arguments of principle 
which would permit us to decide: all gauges, and even the gauge invariant formalism are on 
the same footing; one can work in any gauge. We hope however that in this paper we clearly 
demonstrated that for practical purposes the gauge invariant formalism (or equivalently the 
longitudinal gauge) is by far the most convenient one for treating cosmological perturbations 
and is an insurance against mistakes. 
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